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Abstract. We give a new, simple, dimension-independent definition of the seren- 
dipity finite clement family. The shape functions are the span of all monomials 
which are linear in at least s — r of the variables where s is the degree of the 
monomial or, equivalently, whose superlinear degree (total degree with respect to 
variables entering at least quadratically) is at most r. The degrees of freedom are 
given by moments of degree at most r — 2d on each face of dimension d. We establish 
unisolvence and a geometric decomposition of the space. 



The serendipity family of finite element spaces arc among the most popular fi- 
nite element spaces for parallelepiped meshes in two, and, to a lesser extent, three 
dimensions. For each such mesh and each degree r > 1 they provide a finite ele- 
ment subspace with C'^ continuity which has significantly smaller dimension than the 
more obvious alternative, the tensor product Lagrange element family. However, the 
serendipity elements are rarely studied systematically, particularly in 3-D. Usually 
only the lowest degree examples are discussed, with the pattern for higher degrees 
not evident. In this paper, we give a simple, but apparently new, definition of the 
serendipity elements, by specifying in a dimension-independent fashion the space of 
shape functions and a unisolvent set of degrees of freedom. 

The serendipity finite element space Sr may be viewed as a reduction of the 
space Qr, the tensor product Lagrange finite element space of degree r > 1. The Qr 
elements are certainly the simplest, and for many purposes the best, C° finite elements 
on parallelepipeds. They may be defined by specifying a space of polynomial shape 
functions, Qr{I"')j and a unisolvent set of degrees of freedom, for the unit cube /"■. 
Here / = [—1, 1] is the unit interval, n > 1 is the space dimension, and for / C M"' 
and r > the space Qr{f) is defined as the restriction to / of the functions on 
M"" which arc polynomial of degree at most r in each of the n variables separately, 
so dim Qr{r'') = (r + 1)". In addition to the usual evaluation degrees of freedom 
associated to each vertex, to each face / of /" of some dimension d> 1 are associated 
the degrees of freedom 
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(properly speaking, q ranges over a basis of Qr-2{f), but we shall henceforth dispense 
with this distinction). Since the number of faces of dimension d of /" is 2"^''(|||) , and, 
by the binomial expansion, 

we see that the total number of degrees of freedom coincides with the dimension of 
Qr{I"'). The proof of unisolvence is straightforward, by using the degrees of freedom 
to show, in turn, that u vanishes on the faces of dimension 0, 1, 

The serendipity family in two dimensions is discussed in most finite element 
textbooks, especially the lowest order cases, namely the 4-node, 8-nodc, and 12-node 
rectangular elements. The idea is to maintain the same degrees of freedom as for 
on the boundary of the element, but to remove interior degrees of freedom, and specify 
a correspondingly smaller shape function space for which the smaller set of degrees 
of freedom is unisolvent. In this way, one obtains a space of lower dimension without 
sacrificing C'^ continuity and, hopefully, without much loss of accuracy. For r = 1, 
there are no interior degrees of freedom for Q,., so the space Si is identical to Qi. But 
the serendipity space S2 and S3 have only 8 and 12 degrees of freedom, respectively, 
compared to 9 and 16, respectively, for Q2 and Q3. A possible generahzation to higher 
degree is to keep only the 4r boundary degrees of freedom of Qr and to seek Sr as 
a subspace of Q,. of dimension 4r for which these degrees of freedom are unisolvent. 
This is easily accomplished by taking Sr to be the span of the monomials x*, x''y, |/*, 
and xy', for < i < r + 1, (4r monomials altogether, after accounting for duplicates). 
The resulting finite element space is referred to as the serendipity space in some of 
the literature, e.g., [1]. However, for r > 3 this space does not contain the complete 
polynomial space Vr, and so does not achieve the same degree of approximation as Qr. 
Therefore the shape functions for the serendipity space in two dimensions is usually 
taken to be the span of Vr{I'^) together with the above monomials, or equivalently. 

The degrees of freedom associated to the vertices and other faces of positive codimen- 
sion are taken to be the same as for Qr, and the degrees of freedom in the interior 
of the element can be taken as the moments u i->- Jp uq, q e Vr-4:{I^), resulting in a 
unisolvent set. 

This definition does not generalize in an obvious fashion to three (or more) 
dimensions, but many texts discuss the lowest order cases of serendipity elements in 
three dimensions: the 20-node brick, and possibly the 32-node brick [1, 2, 4, 6], which 
have the same degrees of freedom as Qr{I^) on the boundary, r = 2,3, but none 
in the interior. It is often remarked that the choice of shape function space is not 
obvious, thus motivating the name "serendipity." The pattern to extend these low 
degree cases to higher degree brick elements is not evident and usually not discussed. 
A notable exception is the text of Szabo and Babuska [5], which defines the space 
of serendipity polynomials on the three-dimensional cube for all polynomial degrees, 
although without using the term serendipity and by an approach quite difi^erent from 
that given here. High-degree serendipity elements on bricks have been used in the 
p- version of the finite element method [3] . 
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In this paper we give a simple self-contained dimension-independent definition of 
the serendipity family. For general n, r > 1, we define the polynomial space Sr{I"') and 
the degrees of freedom associated to each face of the n-cube and prove unisolvence. 

2. Shape functions and degrees of freedom 

2.1. Shape functions. We now give, for general dimension n> 1 and general degree 
r > 1, a concise definition of the space 5^(1") of shape functions for the serendipity 
finite element. (As usual, a monomial is said to be linear in some variable Xi if it is 
divisible by but not x'f, and it is said to be superlinear if it is divisible by x?.) 

Definition 2.1. The serendipity space Sr{I^) is the span of all monomials in n 
variables which are linear in at least s — r of the variables where s is the degree of the 
monomial. 

We may express this definition in an alternative form using the notion of the super- 
linear degree. Define the superlinear degree of a monomial p, denoted deg2P, to be 
the total degree of p with respect to variables which enter it superlinearly (so, for 
example, deg2 x'^yz^ = 5) and define the superlinear degree of a general polynomial 
as the maximum of the superlinear degree of its monomials. 

Definition 2.1'. The serendipity space Sr{I"') is the space of all polynomials in n 
variables with superlinear degree at most r. 

It is easy to see that Definitions 2.1 and 2.1' are equivalent, since for a monomial p 
which is linear in / variables, degp = deg2P + I. Surprisingly, neither form of the 
definition seems to appear in the literature. 

Since any monomial is, trivially, linear in at least variables, and no monomial in 
n variables is linear in more than n variables, we have immediately from Definition 2.1 
that Vr{I'^) C Sr{I"') C Vr+n{I^), whcrc Vr{I'^) is the space of polynomial functions 
of degree at most r on J". In fact, Sr{I"') C Vr+n-i{I"'), since the only monomial 
which is linear in all n variables is xiX2 ■ ■ ■ Xn, which is of degree n, not degree r + n. 
In particular, the one-dimensional case is trivial: Sr{I) = Vr{I) = Qr{I), for all 
r > 1. The two-dimensional case is simple as well: Sr{P) is spanned by Vr{I'^) and 
the two monomials X1X2 ^'^d x\x2 (these two coincide when r = 1). Thus we recover 
the usual serendipity shape functions in two dimensions. In three dimensions, Sr{I^) 
is obtained by adding to Vr{I^) the span of certain monomials of degrees r + 1 and 
r -|- 2, namely those of degree r which are linear in at least one of the n variables, 
and those of degree r + 2 which are linear in at least two of them (there are three of 
these X\X'2X^ 5 X^X^X^^ X^X^X'^ except for r — 1, when all three coincide). 

To calculate dimiSr(/") we count the monomials in n variables with superlinear 
degree at most r. For any monomial p in xi, . . . , let J C N„ := {1, . . . , n} be the 
set of indices for which Xi enters p superlinearly and let d > be the cardinality of 
J. Then 

p=(n^|)xgx(na:r), 

where g is a monomial in the d variables indexed by J and each ai, i E J'^ (the 
complement of J), equals either or 1. Note that deg2P — 2d + degq. Thus we 
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Table 1 dim 5^(7") 



may uniquely specify a monomial in n variables with superlinear degree at most r 
by choosing d > 0, choosing a set J consisting of d of the n variables (for which 
there are (^) possibilities), choosing a monomial of degree at most r — 2d in the d 
variables (('^^'') possibilities), and choosing the exponent to be either or 1 for the 
n — d remaining indices (2"^°' possibilities). Thus 

min(n,[r/2j) / \ / _ ^\ 

(2.1) dim 5.(7")= Y: 

Table 1 shows the dimension for small values of n and r. 

2.2. Degrees of freedom. We complete the definition of the serendipity finite ele- 
ments by specifying a set of degrees of freedom and proving that they are unisolvent. 
Let / be a face of 7" of dimension d> 0. Then the degrees of freedom associated to 
/ are given by 

(2.2) [uq, qeVr-2d{f). 

Jf 

Note that Vs{f) is defined to be the space of restrictions to / of Vsil""), so if / is a 
vertex, then Vs{f) — K for all s > 0. In this case the integral is with respect to the 
counting measure, so each vertex is assigned the evaluation degree of freedom. 

In contrast to the degrees of freedom (1.1) of the tensor product Lagrange family, 
for the serendipity family, the degrees of freedom on faces are given in terms of the 
Vsif) rather than Qs{f). For example, for ^6(7^), there arc 6 = d[m'P(i-2x2{l'^) 
degrees of freedom internal to the square. For Sq{P), there arc 6 degrees of freedom 
on each face of the cube and one degree of freedom internal to the cube. Figure 1 
shows degree of freedom diagrams for the spaces 5.(7") for r < 4 and n < 3. 

Theorem 2.2 (Unisolvence). The degrees of freedom specified in (2.2) are unisolvent 
forSrin. 

Proof. First we note that number of degrees of freedom equals the dimension of the 
space. Indeed, there are 2"~'^(^) faces of the cube of dimension d (obtained by fixing 
n — d of the variables to ±1) and, for / of dimension d, dimVr-2d{f) = so the 

total number of degrees of freedom proposed is precisely dim 5.(7") given by (2.1). 
It remains to show that if m G 5.(7") and all the quantities in (2.2) vanish, then 
u vanishes. We do this by induction on n, the case n — 1 being trivial. Let F be 




Fig. 1 Degrees of freedom for <S^(/") for r = 1, 2, 3, 4, n = 1, 2, 3. Dots 
indicate the number of degrees of freedom associated to each face. 



a face of /" of dimension n — 1 and let U = u\f. The f/ is a polynomial in n — 1 
variables with superlinear degree at most r, i.e., U G Sr{F). Moreover, if / is any 
face of F of dimension d, and q e Vr-2d{f), then J^Uq — J^uq — 0. Therefore, by 
the inductive hypothesis, U vanishes. In this way, we see that u vanishes on all its 
faces, i.e., whenever we fix some Xi to ±1. Therefore 

for some polynomial p. Note that deg2 u = 2n + degp. Therefore degp < r — 2n, and 
we make take f = I"' and q = p in (2.2) to find 



Jin 



{l-xl).--{l-xl)p'^0. 



It follows that p = 0. 



□ 



In the proof of unisolvence we established that the degrees of freedom associated 
to a face and its subfaces determine the restriction of u to the face. This is important, 
since it implies that an assembled serendipity finite element function is continuous. 



3. Geometric decomposition 



In this section we give a geometric decomposition of Sr{I'^), by which we mean 
a direct sum decomposition into subspaces associated to the faces. Such a decom- 
position can be used to derive explicit local bases which are useful for the efficient 
implementation of the elements, and also for insight. 

First we introduce some notation. Let Ad{I^) denote the set of faces of the n- 
cube of dimension d and A(/") = 1J^*^q Ad(/") the set of all faces of all dimensions. 
A face / e Ad(/") is determined by the equations xj = cj for j E J where J C N„ is 
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a set of cardinality n — d and each cj e { — 1, 1}. We define the bubble function for / 
as 

the unique (up to constant multiple) nontrivial polynomial of lowest degree vanishing 
on the (n— l)-dimensional faces of /" which do not contain /. Note that hf is strictly 
positive on the relative interior of / and vanishes on all faces of /" which do not 
contain /. For s > we denote by (/") the space of polynomials of degree at most 
s in the d variables Xj, j e J". lid = Q, then is understood to be R. 

To a face / of dimension d we associate the space of polynomials Vf := 'P/_2d(-^")^/- 
By definition, any element q ^ Vf has the form q = pYlj&ji^j where p G 7^/(1"), 
and so deg2 ? < r. Thus Vf C <Sr(/"). The following theorem states that they do 
indeed form a geometric decomposition. 

Theorem 3.1 (Geometric decomposition of the serendipity space). LetSr{I'^) denote 
the serendipity space of degree r, and for each f e Ad{I"'), let Vf — Vl_2d{I^)bf- Then 

(3.1) Sr{n= ^f- 

/eA(7") 

Moreover, the sum is direct. 

Proof Clearly dim^/ = dimP/_2rf(/") = C"/), and so (2.1) imphes that 

dim 5^(7")= dimVf. 
/eA(7") 

Hence, it is sufficient to prove that if p is a monomial with superlinear degree < r, 
then p e Vf. Write p = x" ^ • • • x"" and let 

Jo = { i e N„ I = }, Ji = { i e N„ I = 1 }, J2 = ( Jq U Ji)^ 

denote the sets indexing the variables in which p is constant, linear, and superlinear, 
respectively. By assumption, X]jeJ2 '^i — '^^ expand 

\{l + Xj) + \{l-Xj), J e Jo, 

^7 = {¥^ + '^j) - 1(17 ^j)' 3 e Ji, 

j(l + Xj) + i^(l - Xj) + (1 - x2)r,(x,), J e J2, 

where deg rj = aj — 2, and insert these into p = x"^ • • • x^". We find that p is a linear 
combination of terms of the form 

(3.2) q=llr,{xj){l-x'^) J] 

jeJ' ieJ"uJoUJi 

where J' and J" partition J2. Let / be the face given by 1 =F = 0, j G J" U Jq U Ji, 
where the signs are chosen opposite to those on (3.2). Then dim / = #J' and q = rbf 
with r e 7'/(i""), s = Ejej'(«j - 2) < r - 2dim/. This shows that q e Vf and 
P ^ Z]/ ^ desired. □ 
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